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Abstract. We construct an analytic model for the mass distributions of superclusters 
and clusters in each supercluster. Our model is a modification of the Press-Schechter 
theory, and defines superclusters as the regions that have some overdensity smaller 
than those for usual virialized objects. We compare the mass functions with a catalog 
of superclusters in the Sloan Digital Sky Survey Early Data Release and found that 
they are in reasonable agreement with each other. 
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1. Introduction 

It is now known that our universe is hierarchical. The largest virialized objects in the 
universe are clusters that are made of hundreds of galaxies. Clusters have extensively 
been used to test structure formation models in the universe, through their abundances 
and correlations P3 Ej. 

Superclusters are defined as an ensemble of clusters. They are still assembling, and 
are even larger structure than clusters; superclusters have extents of ~ lOOMpc or even 
larger. Thus superclusters offer us information on late evolution of universe, and also on 
the transition from linear to non-linear regime. Examples of the superclusters include 
A901/902j3i, MS0302+17|11E], C116O4+430 Hj, the Sculptor supercluster^, and the 
Shapley superclusterj^l. Detailed investigations have shown that these superclusters are 
real physical systems where clusters are gravitationally assembling. Detailed studies 
of these superclusters from galaxy distributions and weak lensing also showed that 
the mass distribution of superclusters is in good agreement with the distribution of 
early-type galaxies in superclusters. The physical state of gas in superclusters is still 
unknown, but it is worth mentioning that soft X-ray has been detected in the Sculptor 
supercluster . 

Although most the studies have been limited to those of individual superclusters 
so far, the emergence of wide-field surveys/telescopes, such as the Sloan Digital Sky 
Survey (SDSS) and the Subaru Suprime-cam, begins to make it possible to do statistical 
studies of superclusters. Nevertheless, there is no analytic model of superclusters that 
predicts, e.g., the mass distribution. In this paper, we construct such an analytic model 
of superclusters. Basically, our model is a modification of the Press-Schechter theorv[TT|. 
and predicts not only mass distributions of superclusters but also number distribution 
of clusters in each supercluster. We also compare our model with superclusters in the 
SDSS extracted by Einasto et a/.p2]- Throughout the paper we assume a flat universe 
n^.^ + fi^ = 1, and the Hubble constant h = iJo/(100kms-^Mpc-i) = 0.7. 

2. An Analytic Model 

In this section, we construct the mass function of superclusters and clusters in each 
supercluster using the spherical model[TT]. We regard superclusters as overdense regions 
that are on the course of the spherical collapse. Then superclusters can be defined as 
the regions that have the overdensity Age = Psc/p, where psc is the average density 
within superclusters and p is the mean density in the universe. Although superclusters 
are usually characterized by filamentary structure that is quite different from simple 
spherical overdensity, we consider this model so as to be analytically tractable. The 
assumption, however, may not be so crude given the similar situation that the Press- 
Schechter theory has achieved reasonable success despite its simplicity. 




Figure 1. Sizes of overdense regions (r) and corresponding nonlinear overdensities 

(Age) as a function of dimensionless cosmic time r. Different lines mean different 
values of k: from k — 1 (Dark) to 10 (Light). The vertical dashed line indicates the 
value of T corresponding to a = 1. In this plot, we adopt = 0.3. 



2.1. Spherical Model in Lambda- dominated Universes 

Consider a spherical overdensity with mass M and radius R. Equation of motion of the 
overdensity is simply given by 

^--^ + ^^HoR- (1) 

By integrating this equation, we obtain 

^ 1 fdnV GM 1 

here £^ > for a growing perturbation. We define 

\HqQm J 
_ 2i? 

-no"M-i£o 

T =^/n^Hot, (6) 
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Figure 2. The relation of linear (eq. ^U]) and nonlinear (eq. overdensities at 

a — I. Here we adopt ^Im = 0.3. 



Then equation reduces to 

On the other hand, the cosmological time is related with the scale factor as 

r = ^tu^^/^arccosh (l + 2u;a^) . (9) 

From equation (jS)), we can derive the motion of spherical overdensity as a function of r. 
Figure H shows the evolution of the size of overdense region for various values of n. It is 
shown that the corresponding overdensity Age increases monotonically as time evolves. 
Since the dimensionless time r is related with the scale factor a as equation Q, at fixed 
time a and for fixed Age we can determine the value of k uniquely by solving equations 
(jHl) and Q. Then the linear and nonlinear overdensity is easily derived: 

4c = ^D{a), (10) 
5 

Asc = ^, (11) 

where D{a) is the linear growth rate normalized so that D(a) = a at a ^ 1. From 
this, arbitrary Age is converted to the corresponding linear overdensity 6sc that is an 
essential quantity in computing the mass function. In Figure El we plot an example of 
this relationship. It shows one-to-one correspondence between Age and 6sc ■ 



2.2. Mass Function of Superclusters and Clusters 

To compute the mass function of superclusters, we adopt the Press-Schechter 
formalism |llj. From the linear overdensity 6sc (eq. PHI); the mass function can be 
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Figure 3. Ovcrdcnsities versus total luminosities of superclusters in the SDSS 
EDR 12 . We assumed the average r-band luminosity density pL ~ 2.6 x 
lO^/iL0Mpc-3[IH] to convert into A^c- 



computed as 
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(12) 



where 6sc{z) = 6sc{D{z = 0)/D{z)}, and um is the (hnear) mass variance at 2; = 0. 

The conditional mass function of clusters within superclusters, that are already 
virialized, can be estimated in analog with the extended Press-Schechter formalism |13| 
HU HH] . The number of clusters with mass M that are in a supercluster with mass Mgc 
is given by 
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where Sc{z) ~ 1.68{D{z = 0)/D{z)} is the critical overdensity in the standard Press- 
Schechter formalism. Although the distribution of clusters in each supercluster might 
be affected by physical processes such as global tides and dynamical frictions}^, we 
adopt this form for simplicity. 



3. Comparison with Observations 

In this section, we compare our analytic model with the results of Einasto et a/.^^l who 
made a catalog of superclusters in the SDSS Early Data Release (EDR)[T7j. They 
extracted superclusters from the density field of galaxies with a smoothing length 
10/i~^Mpc, and also clusters from the density field with a smoothing length 0.8h~^Mpc. 
The superclusters and clusters are defined as the regions divided by the threshold density 



An Analytic Model of Superclusters 



6 



(5 = 1.8). Since the SDSS EDR data consists of two slices, the density filed is calculated 
in 2 dimensions only. Thus the total (r-band) luminosity of a supercluster is derived as 

Ltot = -j^Lohs, (14) 

where is the thickness of the slice, D is the average diameter of the supercluster, and 
Lobs is the observed luminosity with the selection function of galaxies corrected. With 
this way 39 superclusters were identified (see Table 3 of [121 )• 

Since the above definition of superclusters is slightly different from theoretical one, 
we need to check what the observed superclusters correspond to in our model. We have 
defined superclusters as the regions that have the overdensity Age, so we compute Age 
for each observed supercluster. Figure El shows A^c for the SDSS EDR superclusters. 
Here we calculated Age as 

PL 

where pl = 2.6 x lO^/iLoMpc"^ is the average r-band luminosity density^HI- This 
coincides with equation (jllj) if the mass-to-light ratio of superclusters is the same as 
that of the universe. As seen in the Figure, the derived value of Age is scattered at 
around 2 (Age = 2.1 ± 0.8). Thus we adopt A^c = 2 as a fiducial value, though we 
also examine the cases with Age 7^ 2 to see how the uncertainty of Age would affect our 
results. 



3.1. Mass Function of Superclusters 

First, we compare the abundance of superclusters in the SDSS EDR field. To do so, 
the volume of the field where superclusters were identified is needed. We simply adopt 
the sum of volumes of North and South slices of the EDR data^2]. However, this 
may underestimate the volume because of the correction outside the slices (eq. [H]). 
To convert mass functions of superclusters to luminosity functions, the mass-to-light 
ratio of superclusters. Again, we assume the fiducial model M/L = AOOHMq/ Lq, and 
estimate the sensitivity by changing M/L from 250hMQ/ Lq to SOO/iMq/Lq. 

Figure |3 plots luminosity functions of superclusters in theory and observations. 
First, the abundance of massive superclusters is sensitive to Age, thus the selection of 
superclusters is quite important in comparing the abundance. Next, the abundance 
of massive superclusters is also very sensitive to ag. Therefore, surveys of massive 
superclusters can be a good probe of ag. Finally, the observed number of superclusters 
shows reasonable agreement with the theory, particularly when we take account of 
several uncertainties including that of the survey volume. At smaller luminosities the 
difference seems large, but this may be ascribed to the incompleteness of observations 
for such small objects. 
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Figure 4. Comparison of luminosity functions of superclusters. The observational 
data {filled triangles with errorbars) are taken from |12| . For the errorbars, only 
statistical errors are included. We adopt the following fiducial model in computing 
theoretical predictions: A^c = 2.0, M/L = AOOhMQ/LQ, = 0.3, and ag = 0.9. For 
each panel, we changed one of these parameters while the other parameters are fixed. 
In all panels, the darker lines mean the larger values. 



3.2. Number of Clusters in Superclusters 

Einasto et al. also estimated the number of clusters in each supercluster. They noted 
that there exists a well-defined lower limit of cluster luminosities that is a function 
of distance to clusters (denoted by rf), particularly for those at larger distances (see 
Figure 6 of P^). Thus we divide the supercluster sample by the distance d. For each 
subsample, we assume the limiting luminosity shown in Figure 6 of [T2j and calculate 
the number of clusters in a given supercluster by integrating equation (fTHj). We also 
convert the observed number of clusters in the EDR field to the total number in the 
similar way as equation (fT^ . 

In Figure we plot the number of clusters within superclusters for different 
distances d. We find that theoretical curves and observations show good agreement 
with each other. It seems that the model overpredicts the number of clusters, but this 
is perhaps because the assumption of the existence of lower limit of cluster luminosities 
becomes inaccurate for lower d (see Figure 6 of Next we examine the dependence 

of several parameters on the number of clusters, which is shown in Figure IHl The 
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Figure 5. Comparison of the numbers of clusters in superclusters. The observational 
data (filled triangles) are taken from |12j . For the errorbars, only statistical errors 
are included. We consider four subsamples of superclusters divided by the distances 
to superclusters; d ^ 200/i^iMpc (ISO/i^^Mpc < d < 250/i~iMpc; upper left), 
d - SOO/i-^Mpc (250/i~^Mpc < d < 350h~^Mpc; upper right), d - 400/i-^Mpc 
(350/i-iMpc < d < 450/i-iMpc; lower left), and d - SOO/i-^Mpc (450/i-iMpc < 
d < 550/i~^Mpc; Zower right). As limiting masses of clusters, we adopt "set 1" in 
Figure 6 of P2|- Different lines mean different values of Age- The other parameters 
are fixed to the fiducial values (see the caption of Figure . 



Figure indicates that the number of clusters in superclusters hardly depends on these 
parameters, though that was already clear from the expression of the conditional mass 
function shown in equation (fTSj) . Thus the model predictions are quite robust, and the 
good agreement means that our model works well. 

4. Summary and Discussion 

In this paper, we have constructed an analytic model of superclusters. Our model 
is based on the Press-Schechter theory, and defines superclusters as overdense regions 
that are on the course of the spherical collapse. Our model allows us to compute the 
mass function of superclusters, and also the number distributions of clusters in each 
supercluster. 

To test the model, we have compared the model with the supercluster catalog in the 
SDSS EDR data. We have shown that the model shows in reasonable agreement with 
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Figure 6. Comparison of the numbers of clusters in superclusters at d 500/i~^Mpc. 
The observational data (filled triangles) are taken from |12j . For the errorbars, only 
statistical errors are included. For each panel, we changed one parameters while the 
other parameters are fixed to the fiducial values (see the caption of Figure . In all 
panels, the darker lines mean the larger values. 



the observation. We have also examined parameter dependences of the mass function 
of superclusters, and found that the mass function of massive superclusters is sensitive 
to CTg. Thus superclusters will be a useful probe of if superclusters are selected more 
rigorously. On the other hand, the number of clusters in superclusters shows only weak 
dependences on parameters. 

Our model would be useful not only in constraining cosmological parameters from 
the abundance of superclusters but also in many applications. One of such applications 
is hot gas in superclusters. Recently it has been claimed that soft X-ray^U] and 
Sunyaev-Zel'dovich effectsjinj in the supercluster regions may be detected, and this 
might suggest the existence of relatively hot gas in superclusters. If this is true, the 
statistical significance of the signals that are caused by the hot gas can be computed 
by using our model, once the we measure the temperature and density structure of 
hot gas. Other application is gravitational lensing; gravitational lensing measures the 
total mass along the line of sight, thus a supercluster outside a cluster may affect, e.g., 
the mass estimate of the cluster. Our model is useful since it allows us to estimate 
the distribution of the host superclusters for a given cluster. In practice, we found 
that the effect is negligible when we consider superclusters with Age ~ 2; assuming 
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a homogeneous sphere, convergence of a supercluster with M = lO^^h ^Mq can be 
estimated as 



where we assumed a source at z ~ 1 and a supercluster at 2 ~ 0.3. However, it 
is of interest to consider the effect more carefully, because the effect is dependent on 
the density profile of superclusters as well as the definition of superclusters. Age- In 
particular, superclusters could have much more significant effect if superclusters have 
centrally concentrated density distributions. 

Although our model looks successful at this moment, definitely we need more tests. 
The more direct way to test the model is to compare it with numerical simulations, since 
we know the initial condition and cosmological parameters in numerical simulations. Of 
course, comparison with observations is also of great importance in order to test our 
understanding of the universe. Now the SDSS data have been becoming public jSOl EI], 
and the quite large survey region of the SDSS would allow us to make three-dimensional 
catalog of superclusters. The follow-up work is currently in progress, and the results 
will be presented elsewhere. 
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